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ABSTRACT 

We present early-time optical through infrared photometry of the bright gamma-ray burst 
GRB 080607, starting only 6 s following the initial trigger in the rest frame. Complemented by 
our previously published spectroscopy, this high-quality photometric dataset allows us to solve for the 
extinction properties of the redshift 3.036 sightline, giving perhaps the most detailed information on 
the ultraviolet continuum absorption properties of any sightline outside our Local Group to date. The 
extinction properties are not adequately modeled by any ordinary extinction template (including the 
average Milky Way, Large Magellanic Cloud, and Small Magellanic Cloud curves), partially because 
the 2175 A feature (while present) is weaker by about a factor of two than when seen under similar cir- 
cumstances locally. However, the spectral energy distribution is exquisitely fitted by the more general 
Fitzpatrick & Massa (1990) parameterization of Local-Group extinction, putting it in the same family 
as some peculiar Milky Way extinction curves. After correcting for this (considerable. Ay = 3.3 ± 0.4 
mag) extinction, GRB 080607 is revealed to have been among the most optically luminous events 
ever observed, comparable to the naked-eye burst GRB 080319B. Its early peak time (trcst < 6 s) 
indicates a high initial Lorentz factor (F > 600), while the extreme luminosity may be explained in 
part by a large circumburst density. Only because of its early high luminosity could the afterglow of 
GRB 080607 be studied in such detail in spite of the large attenuation and great distance, making 
this burst an excellent prototype for the understanding of other highly obscured extragalactic objects, 
and of the class of "dark" GRBs in particular. 

Subject headings: Gamma-ray burst: individual: 080607 — dust, extinction 



1. INTRODUCTION 

The most extreme gamma-ray bursts (GRBs) have of- 
ten been the most illuminating — both literally and 
figuratively. The enormous isotropic-equivalent en- 
ergy of GR B 971214 (redshift z = 3.43, ^isn = 
3 X lO^'^ erg: iRamaprakash et al.lll998t lOdewahn et al.l 



Il998t iKullarni et all "9981 ) emphatically demonstrated 
the need for collimation to bring the energy budget of 
long-duration GRBs within physically reasonable values. 
Observations of the mag 9 optical fl ash of GRB 990123 
(z = L61, ^isn = 3.4 X lO^'' erg; lAkerlof et all [l99l 
iKulkarni et al.lfT999( ) anticipated the utility of GRBs to 
probe the high-redshift universe: similar events would be 
easily detectable even at z > 6. This possibility was first 
vindi cated by GRB 050 9 04 (z = &.2<^ E^a = 1.2 x 10^'* 
erg; iKawai et all l2006t ISugita et all l2008[ ). which for 
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three years remained the most distant GRB known and, 
at the time, was also the most luminous optica l tran- 
sient observed in the Universe (jKann etaD|2Q03). The 
latter record has since been surpassed dramatically by 
GRB 080319B (z = 0.937, E;^^ = 1-3 x lO ^^ erg), whose 
optical afterglow p eaked at F » 5 mag (jRacusin et al.l 
120081 iBloom et al.ll2009l : IWozniak et al.ll2009t ). The cur- 
rent record for the bolometric isotropic-equivalent en- 
ergy is held b y the J^ermz burst GRB 080916C (Ej^ p = 
6.5 X lO^-* erg; lAbdo et al.ll2009t iGreiner e t al 2009). 

Joi ning this list of re cord setters is GRB 080607 (z — 
3.036: iProchaska et al.ll2009D . with E-.^^ = 1.87 x 10^^ erg 
(|Golenetskii et al.ll2008[ r" This event is remarkable not 
only for its intrinsic properties, but also because of its 
unusual environment: a Keck spectrum obtained start- 
ing only 20 min after the burst (Prochaska et ah, 2009[) 
reveals that the sightline penetrates a giant molecular 
cloud in the host galaxy, obscuring the rest-frame visi- 
ble light by Av ~ 3 mag of extinction (or ~ 6 mag at 
1600 A, corresponding to the observed R band) before it 
even began its journey through intergalactic space. In 
spite of this extreme attenuation, the event was bright 
enough to be detected by small optical telescopes for over 
an hour. 

The spectroscopic pro perties of this ey e nt ha ve been 
previously discussed bv IProchaska et al] ()2009() . along 
with a preliminary analysis of its extinction properties; 
further analysis of the spectra was also presented by 

Dust extinctio n is limited or absent for the vast majority of 
well-studied GRBs HSchadv et al.l 2007; Kann et al. 201fl,), and no 
other GRB displays firm evidence for molecular lines. 
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ISheffer et al.l ()2009l ). In this paper, we analyze several 
other aspects of this burst, from the prompt emission 
(and simultaneous optical detection) through a late-time 
search for the host galaxy, and we present a significantly 
expanded discussion of its extinction properties. In Sj2] 
we describe our early-time multicolor observations of the 
afterglow with several different robotic telescopes. We 
analyze the optical light curve in j|3.1H3.2l and show no 
correlation between the prompt emission behavior and 
the early optical observations, starting at only 6 s post- 
trigger in the host frame, and we present limits on color 
variations at early times. In §3.3H3.4l we examine in more 
detail the combined photometric and spectroscopic spec- 
tral energy distribution (SED) and place our final con- 
straints on the host-galaxy extinction properties, demon- 
strating the firm detection of a 2175 A bump, the highest- 
redshift detection of this signature to date. The X-ray 
light curve is analyzed in i )3.5l to search for evidence of 
dust scattering in the host at these wavelengths. In §4.11 
we place GRB 080607 and its environment in the context 
of other GRBs, both ultraluminous and bright events like 
GRB 080319B as weh as the poorly understoo d class of 
extremely dark bursts (I Jakobsson et al.l l2004t ) . Finally, 
in m.2\ we attempt to explain the origin of the burst's 
extreme luminosity. 

2. OBSERVATIONS 
2.1. Swift 

GRB 080607 triggere d the Burst Alert Telescope 
(B AT, Barthclmy ct al] l2005l ) on the Swift satellite 
(jGehrels et al.l 120041) at 06:07:27 on 2010 June 6 (UT 
dates are used throughout this paper; times are refer- 
enced to this trigger time, although it is important to 
note that there was significant emission before this trig- 
ger). The light curve (Figure [T]) is spiky and erratic, 
exhibiting a dominant peak at ~ 4 s as well as numerous 
other, fainter peaks ranging from a few seconds before the 
trigger out to ~130 s after, when the signal falls below 
the background level. Swift slewed immediately to the 
source and begai i pointed observat ions with the X-ray 
Telescope (XRT. iBurrows et al.l ^051 at 82 s, followed 
by obser vations with the U ltraviolet Optical Telescope 
(UVOT; iRoming et aIll2005D beginning at 100 s. Obser- 
vations continued until 1049 s, after which Swift slewed 
away temporarily, returning to the field at 4226 s. From 
then, observations continued intermittently over the next 
four days, after which the X-ray flux was too faint for 
Swift to detect. 

The BAT observations were processed using the Swift 
HEAsoft 6.5 software package via the burst pipeline 
script, batgrbproduct. We calculated spectral parame- 
ters both direc tly and u s ing th e Bayesian formalism de- 
scribed by B utler et al.l (l2007l ). Fit to this burst alone, 
a GRBM ((Band et al.lll993l ) model provides no signifi- 
cant improvement over a basic power-law fit over BAT's 
15-350 keV energy range (photon index F = 1.16), 
suggesting a peak energy above the BAT range. Us- 
ing the Bayesian estimate of £^p eak,obs = 902j^4^p° keV 
and the measured redshift 3.036 (jProchaska et al.|[2009[ ). 
we estimate a broad-band isotropic-equivalent energy of 
Eiso = 2.8tli >^ lO'^'^ c^rg. These values place GRB 080607 
second in Eiso rank among all Swift GRBs to date and 
in the same regime as extreme events as GRBs 080319B 
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Fig. 1. — Gamma-ray light curve of GRB 080607 (at a combi- 
nation of 128 ms and 1 s binning), sliowing the bright initial pulse 
complex followed by an additional series of pulses lasting for the 
next several minutes. Original data for the light curve are taken 
from the Swift Burst Analyser IIEvans et al.|[2010l') . 

and 990123. 

GRB 080607 was observed by other satellites as well 
(Konus-Wind and Super- AGILE), enabling a precise 
measurement of the spectral parameters. An in-depth 
analysis of the Konus data will be presented in future 
work by Sbarufatti et al., bu t preliminary calculations 
from iGolenetskii et all |2008| give the following values: 
i^peak.obs = 394t^8 keV and E-,,, = 1.87t°;}J x lO^^ erg. 
These are at the low end of, but generally consistent with, 
the Bayesian Swift result, and confirm that GRB 080607 
was among the most luminous and intrinsically hardest 
(highest i?pcak,rcst) GRBs observed by any satellite. 

The X-ray afterglow was detected throughout the ob- 
ser vations; XRT data were reduced by the procedures 
of iButler fc Kocevskil (l2007h . The UVOT afterglow, 
by contrast, is only marginally detected i n the earliest 
epoch , and only in White and V filters (jSchadv et al.l 
l2008f l. Both filters are heavily impacted by damped 
Lyman-a absorption at z = 3.036, and so are not used 
in our analysis. 

2.2. ROTSE 

The ROTSE-III (Robotic Optical Transient Search Ex- 
periment) array is a worldwide network of 0.45 m robotic, 
automat ed telescopes, built for fast responses to GRB 
triggers (jAkerlof et al.ll2003[) . ROTSE-IIIb, located at 
the McDonald Observatory, Texas, responded immedi- 
ately to the i nitial Gamma-ray Bu rst Coordinate Net- 
work (GCN, [Barthelmv et al.lll995( ) alert. The first im- 
age started at 06:07:49.0 UT, 22.0 s after the burst, 
clearly detecting a bright afterglow at the XRT position 
in this exposure. All ROTSE-III images were processed 
with our custom RPHOT photo metry program based on 
the DAOPHO T (lStetsonlll987D PSF-fitting photometry 
package (Ouim bv et al.ll2006b[ ). The unfiltered thinned 
ROTSE-III CCD has a peak sensitivity in the wavelength 
range of R band. The ROTSE magnitudes were thus ad- 
justed using the median offset from the USNO Bl.O R 
band measurements of selected field stars. Observations 
are presented (along with photometry from all other tele- 
scopes, below) in Table [T] 
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2.3. Super-LOTIS 

Super-LOTIS (Livermore Optical Transient Imaging 
System) is a robotic 0.6 m telescope dedicated to the 
search for o ptical counterparts of GRBs (.Williams et al.l 
12004 I2Q08D . The telescope is housed in a roll-off-roof 
facility at the Steward Observatory Kitt Peak site near 
Tucson, AZ. Super-LOTIS triggered on GRB 080607 and 
began observations at 06:08:03 (35 s after the trigger), 
acquiring a series of frames in the R band. The images 
were reduced and photometry performed using standard 
techniques, calibrated relative to nearby Sloan Digital 
Sky Survey (SDSS) standard stars. 

2.4. KAIT 

The Katzman Aut omatic Imagin g Telescope (KAIT) 
at Lick Observatory (" Li et al.ll2003l ) also responded au- 
tomatically to the Swift alert and began taking observa- 
tions, the first starting at 06:09:25, 118 s after the BAT 
trigger. The KAIT filter sequence consists of a series 
of unfiltered observations, followed by a cycle through 
/, and unfiltered exposures. The optical afterglow 
was detected in all filters, although it is quite faint in 
the V band. Following this sequence, a series of un- 
filtered and /-band exposures was manually added, al- 
though the afterglow was not detected in the / band and 
only marginally detected in our unfiltered exposures at 
that time (even after stacking). 

Images were reduced using standard techniques. This 
left a small amount of residual on the background sky, 
which was removed by subtraction of an illumination 
frame. We used aperture photometry to measure the 
afterglow flux, calibrating relative to SDSS stars in the 
field transformed to the Johnson/Cousins system using 
the equations of Lupton (2006). Th e clear-band e xpo- 
sures were calibrated to the R band ()Li et al.ll2003[ ). 

2.5. PAIRITEL 

The robotic Peter s Automatic Infrar ed Imaging Tele- 
scope (PAIRITEL; iBloom et all 1200 6^) consists of the 
1.3 m Peters Telescope at Mt. Hopkins, AZ — for- 
merly used for the Tw o Micron All Sky Survey (2MASS; 
iSkrutskie et al.ir2006f ) — refurbished with the southern 
2MASS camera. PAIRITEL uses two dichroics to image 
in the infrared (IR) J, H, and Kg filters simultaneously 
every 7.8 s. 

PAIRITEL responded to the initial BAT alert and 
slewed immediately to the source. Observations began 
at 06:08:44, 77 s after the trigger, and continued for the 
next 1.3 hr until the source reached its hour- angle limit. 

The early-time (<0.3 hr) raw data files were processed 
using standard IR reduction methods via PAIRITEL 
Pipeline III (Klein et al., i n prep) and resampled using 
SWarp (|Bertin et al.l |20Q2[ ) to create final 1.0" pixel"^ 
images for final photometry. Due to changing sky condi- 
tions that complicated the otherwise superior Pipeline III 
reductions in the Kg band as the source approached the 
horizon, the remainder of the raw data were reduced us- 
ing an older pipeline which utilized a "dark bank" which 
more robustly handles flat-fielding in such cases. 

PAIRITEL 's standard observing cycle is to take three 
7.8 s exposures in immediate succession at each dither 
position. While the early afterglow is detected in 
even the shortest 7.8 s frames, for signal-to-noise ratio 



(S/N) and calibration considerations, we report 23.4-s 
"triplestacks" (the median of all three images at each 
dither position) as our shortest exposures. These images 
were further binned at successively later times to further 
improve the S/N. 

Aperture photometry was performed using custom 
Python software, utilizin g Source Extractor (SExtractor; 
iBertin &: Arnoutslll996D as a back end. Four calibration 
stars present in all images were chosen based on bright- 
ness, proximity of nearby contaminating sources, and lo- 
cation relative to bad pixels. The optimal aperture of 
~3" radius was determined by minimizing the absolute 
error relative to 2MASS magnitudes of our four calibra- 
tion stars. 

Calibration was performed by redetermining the zero- 
point for each image individually by comparison to 
2MASS magnitudes using these four stars. The result- 
ing statistical uncertainty in the zero-point is negligible 
relative to other sources of error. Additiona l, systematic 
source s of error are addressed in detail by iPerlev et al.l 
lj2010f ): we use a similar procedure here to determine the 
total uncertainty of each point. 

2.6. P60 

The robotic P alomar 60 inch telescope (P60; 
ICenko et al.l I2006D automatically responded to the 
Swift trigger for GRB 080607, executing a predefined 
sequence of observations in the Kron R and Sloan i' 
and z' filters beginning 174 s after the burst trigger 
time. Individual images were reduced in real time using 
standard IRAF^^ routines. The images were calibrated 
with respect to several dozen field star s from the SDSS 
Data Release 7 (jAbazajian et al."2009'), using the filter 
transformations of Jordi et al.. (2006.) for the Kron R 
fiher. 

2.7. Keck Spectroscopy 

We initiated spectroscopic observations of the after- 
glow w ith the Low Re solution Imaging Spectrometer 
(LRIS; lOke et al.llT995[) on the Keck I 10 m telescope 
at 13 min after the Swift trigger, although due to poor 
guiding this first frame was not usable. The first ex- 
posure used in our analysis began at 20.1 min follow- 
ing the trigger. Several additional exposures were taken 
over over the next two hours using the B600 grism and 
both the R400 and R1200 gratings; our final observations 
span a wavelength range of 3000-9000 A. Observations 
were fiux-calibrated relative to the spectroscopic stan- 
dard HZ 44. More details on these spe ctroscopic obser- 
vation s and our reductions are given by iProchaska et al] 
(|2009l) . 

2.8. Keck Host-Galaxy Imaging 

The field of GRB 080607 was imaged in several deep 
integrations at Keck in various optical/IR filters from 
g through Kg- None of these integrations resulted in a 
secure detection of the host galaxy, with the possible ex- 
ception of the g band, in which a faint source is detected 
at about 3cr above background at the afterglow position 

IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association for Research in 
Astronomy, Inc., under cooperative agreement with the National 
Science Foundation (NSF). 
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TABLE 1 
Photometry of GRB 080607 



Telescope/GCN 




Filter 


Exp. time 


Mag.'' 


Flux'^ 




sec 




sec 




MJy 


PAIRITEL 


89.0 


J 


23.4 


13.766 ± 0.107 


5048.6 ± 475.4 


PAIRITEL 


89.0 


H 


23.4 


12.050 ± 0.109 


15657.3 ± 1501.3 


PAIRITEL 


89.0 


Ks 


23.4 


10.750 ± 0.139 


33681.7 ±4058.1 


KAIT 


188.0 


clear 


20.0 


17.501 ± 0.055 


363.9 ± 18.0 


KAIT 


158.0 


I 


20.0 


16.582 ± 0.094 


635.8 ±52.7 


KAIT 


128.0 


V 


20.0 


17.538 ± 0.142 


437.9 ±53.7 


ROTSE 


24.5 


clear 


5.0 


14.920 ± 0.040 


3920.9 ± 141.8 


SuperLOTIS 


40.8 


R 




15.060 ± 0.059 


3446.5 ± 182.8 


P60 


234.8 


R 


60.0 


17.524 ±0.014 


356.3 ±4.6 


P60 


406.1 


i 


60.0 


18.359 ± 0.022 


186.4 ±3.7 


P60 


491.8 


z 


60.0 


18.694 ± 0.092 


135.6 ± 11.0 



Note. — Contains only the first data point in each filter taken by each 
telescope. A complete table of photometry will be published online. 
^ Exposure mid-time, measured from the Swift trigger (UT 04:08:54). 
^ Observed value, not corrected for Galactic extinction. 
Corrected for Calactic extinction {Eb-v = 0.07 mag). 



TABLE 2 
Host-Galaxy Limits 



Instrument 


Obs. date 
(UT) 


Exp. time 
(s) 


Filter 


5(T limit 
(mag) 


Keck I / LRIS 
Keck I / LRIS 
Keck I / NIRC 


2009-02-19 
2009-02-19 
2009-05-31 


2490 
2220 
3600 


9 
I 

Ks 


> 27.3 

> 25.3 

> 21.6 


Note. — 5cr 


limiting ma 


gnitudes on 


a host 1 


galaxy at 



the afterglow position from our ground-based optical and 
IR observations at the Keck Observatory. 



{g = 27 A ± 0.3 mag). A log of our ground-based host 
observations is reported in Table [21 

The host galaxy is, however, well detected at 1.6 fim in 
a deep Hubble Space Telescope (HST) image using WF3, 
as well as in both of the warm Spitzer IRAC channels 
(3.6 /xm and 4.5 /xm). The extreme optical faintness of 
this system, while partially due simply to its high red- 
shift [z = 3.036), makes this galaxy of particular interest: 
determination of the redshift would be exceptionally dif- 
ficult using traditional field-survey techniques, illustrat- 
ing the unique ability of GRBs to select and study opti- 
cally faint galaxies at high redshift. Further discussion of 
the host galaxy, including detailed analysis of both the 
ground- and space-based imaging, will be presented in 
upcoming work by Chen et al. 

3. ANALYSIS 
3.1. Optical Light Curve 

The multi-band light curve of GRB 080607 is plotted 
in Figured! After an initially slow decay the light curve 
steepens (decay index a — 1.65, using the convention 
F cx before flattening out at 1000 s to a temporar- 
ily flat decay. This slow decay lasts for approximately 
another hour before fading again, becoming undetected 
in our final KAIT exposures. 

The lig ht curve was fitted using the techniques de- 
scribed bv lPerlev et al.l (|2010[ ) and previous works by our 
group, modeling the light curve as the sum of several bro- 
ken power laws. Our temporal coverage of this event is 
limited (en ding at 10^ s), mak i ng the analysis simple: we 
employ two iBeuermann et al.l ()1999l ) broken power laws. 



one to describe the early behavior and a second to de- 
scribe the later flattening. Because we do not detect the 
rising phase of the afterglow, the pre-break index of the 
first power-law component is not usefully constrained by 
our data and is fixed arbitrarily to —0.5; the falling phase 
of the second (late) component is similarly not well con- 
strained and we fix the post-break decay index to 2.0, the 
shallowest value required for consistency with our upper 
limits in this model. 

Modest but significant color change has been pre- 
viously observed in early-time GRB afterglows (see 
iPerlev et af] |20Q8| or iBloom et all |2009| for two promi- 
nent examples), a possibility which we model by allow- 
ing the intrinsic spectral power-law index f3 {F (x v^^) 
to vary between components or across breaks. However, 
in the case of GRB 080607, any such color change is 
not significant: the change in intrinsic index between the 
fast-decay and fiat components is only A/3 = 0.07 ± 0.07 
and only modestly improves the goodness of fit. There- 
fore, for simplicity we assume no color change during our 
observations of this burst. 

3.2. Absence of Optical/ High- Energy Correlations 

Our optical follow-up observations of this burst begin 
extremely early. The ROTSE coverage begins at only 
21 s after the BAT trigger, corresponding to less than 6 s 
in the GRB rest frame. The prompt emission was still 
extremely active at this time: at least five major gamma- 
ray fiares occurred during our optical observations, the 
last of which was also caught at X-ray wavelengths by 
the XRT. PAIRITEL, KAIT, and SuperLOTIS were aU 
observing during this last flare. 

Even in this rich overlapping dataset, there is no cor- 
relation visible between the optical and high-energ y light 
curves of the type seen by , for example, [V cstran det al.l 
(|200l . iBlake et all (l2005l) . and lBeskin et al.l (l2010l) . In 
Figure [3] we overplot the gamma-ray. X-ray, and optical 
light curves using the same relative scaling. In spite of 
the erratic, flaring high-energy behavior, we see no sign of 
significant deviation of the optical light curves from their 
smooth power-law behavior at any point. Th is is consis- 
tent with other ROTSE- foUowed bursts fe.g. jYost et al.l 
l2007t iRvkoff et"alll2009f ) and provides another clear ex- 
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Fig. 2. — Multi-band light curve of GRB 080607 from a variety of ground-based telescopes as well as the SwiftXRT, fit to a sum of 
two broken power laws. (XRT data are fit to a sum of two unbroken power laws.) Magnitudes are in the Vega system and (with the 
exception of R and Ks) have been shifted as indicated for clarity; these magnitudes are not corrected for Galactic extinction (which is 
nearly insignificant) or host extinction (which is very large). The afterglow initially fades slowly, then steepens; it briefly levels out at 10^ s 
before breaking again and is not detected after 5000 s. The late-time ij-band limits are from Rumyantscv & Pozancnko (200<:|). We use 
the BAT trigger time for tg, which corresponds to the start of the largest prompt-emission pulse; using the start of gamma-ray emission 
instead does not significantly change the qualitative results. 

ample of a burst whose afterglow behavior is clearly di- 
vorced from that of the prompt emission. 

The lack of even modest influence of the prompt emis- 
sion on the afterglow may initially seem surprising: even 
if truly prompt (internal-shock) emission is absent in this 
band, one might expect that some of the energy being 
released so liberally by the central engine might end up 
in the external shock, causing a less dramatic but still 
observ able rebrightening of t he afterglow (a refreshed 
shock; iPanaitescu et al.lll998[ ). We note, however, that 
despite the intense flaring shown in Figure [31 this emis- 
sion is actually dwarfed by an earlier episode: the initial 
pulse of the prompt emission (see Figure [T]) exceeds any 
of the later spikes by an order of magnitude in both inten- 
sity and energy, and it is this initial pulse that dominates 
the energetics of the burst. The later flares are much 
more modest by comparison, so even presuming direct 
input from outflow (revealed by the prompt emission) to 
external shock (revealed by the afterglow), the absence 
of further brightening is not necessarily surprising. 

3.3. Spectral Energy Distribution 

Our light-curve fits naturally provide values for the af- 
terglow flux in each filter at any given time, allowing us 
to model the SED at any time during our observations. 
Because of the absence of significant color change, the 
choice of extraction epoch is arbitrary; 300 s is chosen in 
this case (when all ground-based telescopes were observ- 




60 100 
t(s) 

Fig. 3. — Early-time optical and gamma-ray/X-ray light curve of 
GRB 080607, demonstrating the extremely early peak Ucst S; 6 s 
as well as the absence of any visible correlation between the optical 
and high-energy light curves. (Optical fluxes in other bands have 
been scaled to match the R band; the gamma-ray light curve is 
scaled to match the X-ray curve.) Symbols are the same as in 
Figure [2] 



ing and the afterglow was still bright enough to be well 
detected in all bands). 
The Keck optical spectroscopy covers a wide range of 
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wavelengths and was carefully flux-calibrated: photomet- 
ric standard stars were observed immediately after our 
observations at similar airmasses and the night was pho- 
tometric throughout. Accordingly, we couple our spec- 
trum to the photometry to improve the precision of our 
broadband modeling. 

The optical spectrum is replete with lines from a vari- 
ety of elements and molecules at the host-galaxy redshift 
of 3.036. The analysis of these l ine fe ature s is discussed 
extensively by iProchaska et al.l ()2009D and ISheffer et al.l 
(p009D, and we will not repeat it here; our primary in- 
terest is in the continuum. Although the contribution of 
absorption lines is usually ignored in GRB photometric 
dust modeling, the lines in the spectrum of GRB 080607 
are so abundant and so strong that ignoring them would 
create systematic errors significantly larger than our pho- 
tometric uncertainties in both the spectrum itself and in 
the broadband photometry. In addition, nearly the en- 
tire spectrum at wavelengths shorter than ^6900 A is 
affected by a forest of weak lines from vibrationally ex- 
cited further complicating the analysis. 

Fortunately, we are able to correct for these effects. We 
use the line list presented in Table 1 of IProchaska et aD 
(|2009) to identify all regions of the spectrum affected 
by ionic lines, including the entire spectrum blueward of 
5400 A, which is affected by the host damped Lyman- 
a and the Lyman-a forest. In addition, the spectrum is 
corrected for the subtler but more w idespread a,bsorp - 
tion using the model developed bv ISheffer et al.l (|2009| ). 
We then fit a sixth-order polynomial to the ionic line- free 
regions of this corrected spectrum to create a continuum 
model and perform synthetic photometry using both the 
model spectrum and the observed, uncorrected spectrum 
(and take the ratio) to calculate an adjustment factor 
with which to convert the observed (line-affected) fluxes 
to continuum (line-free) fluxes for each of our broadband 
filters covering the optical spectrum (_R, /, i, and V] 
we assume the line contribution is small further to the 
red). We also wish to use the flux-calibrated spectrum 
itself in later analysis, so we scale the spectrum to the 
photometric SED extraction epoch of 300 s (the scale 
factor is determined by the value that minimizes foi' 
our extinction fits; see §3.3. 1|) and bin the flux in blocks 
of 200 A (excluding line-affected regions). Uncertainties 
are determined by combining the statistical uncertain- 
ties from the spectrum with a systematic term of 3% per 
bin to incorporate any uncertainty in the flux calibra- 
tion (10% is used for <5500 A and >9000 A, which are 
especially uncertain.) Using this technique, we generate 
a line-corrected narrowband SED spanning 5400-9200 A 
to complement our line-corrected photometry. The af- 
terglow fluxes from direct and synthetic photometry are 
presented in Tables [3] and HI respectively. 

3.3.1. Extinction Fitting 

The combined photometric and spectroscopic SED is 
plotted in Figure |4l It is immediately evident that this 
curve is unlike almost any other GRB SED that has been 
observed in detail to date. First, the color is extremely 
red: a power-law fit to the broadband photometry would 
give a spectral slope (F cx z/"^) of /3 « 3, at odds with 
the theoretically expe cted value of = 0.5-1.2 for an 
early fading afterglow (jSari et al.lll998D . Second, it is not 
monotonic: the flux drops sharply from the Kg band until 



V rest (A) 

6000 4000 2000 1000 



14 


1 1 ' 1 ' 1 ' 

- H 






16 








18 




R 


FM ■ ^ 


1 III 'J r'isb%: 




CCM '■■1 


1 ij i 


l\ 

' "-.\\ 


20 


Milky Way 

- - LMC 
SMC 




22 


- CoiTected for Galactic Ay = 0.07 






Host redshift z = 3.036 







20000 10000 6000 

>^eff(A) 



Fig. 4. — The combined photometric and spectroscopic SED of 
GRB 080607 fitted with several different extinction models. (FM 
= Fitzpatrick & Massa 1990, CCM = CardeUi, Clayton, & Mathis 
1989). Black error bars indicate broadband photometry; blue er- 
ror bars show the binned pseudo-photometry as derived from the 
Keck spectrum and corrected for line absorption (including 
The light-grey line shows the spectrum (mostly unbinned and in- 
cluding all lines). Several different extinction fits are shown; only 
the general FM model (solid black) is an acceptable fit to the data. 
The SMC curve shown is a fit to the IR data only (an SMC fit to 
all data converges to Ay = mag). 

~2200 A in the rest frame before actually recovering, 
showing a local maximum at ^^1600 A before falling again 
further to the blue. 

These properties are immediately recognizable as sig- 
natures of dust extinction, and particularly of Milky- 
Way like extinction with its broad 2175 A absorption 
band. This strong extinction imprint, in combination 
with our high-S/N afterglow observations spanning the 
entire optical/near-IR window, permits analysis of the 
rest-frame UV extinction properties at a level of detail 
that is almost never possible with GRBs (or indeed, with 
any other technique at this redshift range). 

To constrain the dust properties, we initially followed 
the standard procedure for GRB extinction measure- 
ments by fitting the average Milky Way (MW), Large 
Magellanic Cloud (LMC), and Small Magellanic Cloud 
(SMC) curves, assuming an intrinsic power-law sectrum. 
(Here, and elsewhere unless otherwise specified, the in- 
trinsic spectral slope over the optical range is fixed at 
/3 = 0.7. Fortunately, because the amount of extinc- 
tion for this burst is so large, deviations from this as- 
sumption do not significantly affect our results, except to 
slightly increase the uncertainties in the derived parame- 
ters , as we will discuss in ii'SA\ ) In all three cases we use 
the iFitzpatrickl (|1999t ) parameterization of Local-Group 
extinction as implemented in the GSFC IDL package, 
with Ry fixed to their average value for each galaxy; for 
SMC extinction we use the Fitzpatrick parameters from 
iGordon et all (l2Q0l (SMC bar average). SMC extinc- 
tion is ruled out (it converges to Ay — mag with x^/dof 
= 1159/24), as it rises steeply to the far-UV (FUV) and 
does not allow for the 2175 A bump feature that is so 
prominent in our data. The LMC and MW curves fit 
the data much better, but nevertheless they are not sta- 
tistically acceptable either. Both curves are too fiat in 
the observed IR; the MW curve also significantly overes- 
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TABLE 3 

Model Fluxes at t = 300 s 





^obs 


w . 

^ obs 




AA,Gal 




'^cont 


^ COllt 


^A.host 




(A) 




(mag) 


(mag) 


(mag) 


(AB mag) 


(MJy) 


(mag) 


X-ray 


12.4 


34.04 














V 


5505 


94.25 ± 11.4 


18.98 ± 0.13 


0.07 


0.60 


18.29 ± 0.13 


174.6 ± 21.1 


6.23 


R 


6588 


200.6 ± 12.6 


17.97 ± 0.07 


0.06 


0.21 


17.88 ± 0.07 


256.6 ± 16.1 


5.76 


i 


7706 


260.5 ± 16.6 


17.86 ± 0.07 


0.05 


0.06 


17.75 ± 0.07 


287.5 ± 18.3 


5.88 


I 


8060 


214.0 ±19.1 


17.64 ± 0.10 


0.04 


0.07 


17.97 ± 0.10 


236.0 ± 21.0 


6.05 


z 


9222 


242.5 ±21.4 


17.96 ± 0.10 


0.03 





17.91 ± 0.10 


249.8 ± 22.0 


6.12 


J 


12350 


867.4 ±51.4 


15.66 ± 0.06 


0.02 





16.54 ± 0.06 


883.1 ± 52.3 


4.94 


H 


16620 


2296 ± 131 


14.12 ± 0.06 


0.01 





15.48 ± 0.06 


2322 ± 132 


4.20 


Ks 


21590 


5866 ± 337 


12.64 ± 0.06 


0.01 





14.47 ± 0.06 


5911 ± 339 


3.30 



Note. — Broadband afterglow fluxes as determined by the light-curve model, interpolated to 
t — 300 s after the trigger. Observed magnitudes are not corrected for Galactic extinction and are 
in the Vega system (except for the SDSS i and z filters, which are given in the SDSS filter system.) 
Continuum magnitudes and fluxes have been corrected for both Galactic extinction (from NED) 
and line absorption (calculated using our optical Keck spectroscopy). 



TABLE 4 

Binned, 
Line-Interpolated 
Keck Spectroscopy 



(A) (/.Jy) 



5448.56 128.24 ± 0.94 

5670.14 146.45 ± 1.79 

5842.98 160.54 ± 2.48 

6112.85 201.91 ± 2.52 

6235.50 196.32 ± 1.74 
6476.67 211.72 ± 2.45 
6776.25 233.69 ± 1.48 
7099.05 240.99 ± 1.72 
7281.72 238.65 ± 1.04 

7483.82 231.69 ± 0.97 
7774.03 213.43 ± 2.41 
7997.30 195.15 ± 0.84 

8181.51 187.27 ± 0.85 
8411.10 180.74 ± 2.26 
8607.66 165.18 ± 1.02 
8759.01 153.90 ± 6.60 

8892.83 151.97 ± 4.61 
9044.91 162.99 ± 2.90 
9175.62 158.64 ± 2.36 



Note. — Uncer- 
tainties are photomet- 
ric only and do not 
include any systematic 
term. Fluxes are cor- 
rected for H2 absorp- 
tion in the host galaxy 
and for Galactic extinc- 
tion. 

timates the strength of the 2175 A bump. 

This should not be a surprise: even within our own 
Galaxy a significant diversity of extinction laws is ev- 
ident. The majority of observed Galactic sightlines are 
consistent with variation in a single parameter Ry , which 
describes the relative "greyness" (wavelength indepen- 
dence) o f the extinction at optical through UV wave- 
lengths (ICardeUi et alj 119891 : hereafter CCM). A small 
number of sightlines in the MW (and all sightlines within 
the LMC and SMC) require additional parameters to fit 
accurately. A more general Local-Group extinction law, 
developed by Fitzpatrick & Massa (1990; hereafter FM), 



is able to fit essentially all local sightlines by adding an 
additional family of parameters: C2 for further variations 
in steepness in the UV, C3 for the strength of the 2175 A 
bump, 7 for the bump's width, and C4 for the strength of 
the FUV rise. (The parameter Ci is also present in prin- 
ciple, but it is essentially degenerate with C2 and Ry, and 
in practice it is fixed based on those values. In addition, 
the parameter xq describes the central wavelength of the 
2175 A bump, but it has not been conclusively shown to 
vary and is fixed to the average value.) 

We first attempted to fit using the general FM law 
(joined to the standard CCM law in the rest-frame op- 
tical with a spline), leaving all parameters free (except 
Ci and Xq as described above). Unfortunately, because 
our observations do not extend far enough into the rest- 
frame optical to properly constrain the optical/IR ex- 
tinction properties independent of the UV, the Ry pa- 
rameter is effectively unconstrained in this case. Fortu- 
nately, Rv and C2 also are tightly correlated locally and 
can be t ied together — us ing, for example, t he corre- 
lation of Fitzpatrick (1999) (linear) or that of 'Reichart| 
([2001) (quadratic, allowing for the optically flat, steep- 
UV SMC-like curve). Both correlations give acceptable 
(and very similar) fits to our data, and the Fitzpatrick- 
constrained curve is shown in Figures [Hini 

We also attempted a range of non-FM models, such 
as t hose of Calzett j e^al. (2000), Maiolin o et al.. (2004), 
and iGaskelfet all §0M)- These curves all lack the 
2175 A bump and do not fit the data. In addition, we 
tried to fit th e multi-parametric extinction curve from 
iLi et all (|2008'), which can incorporate the 2175 A bump 
and gives a fairly reasonable fit (however, the ci parame- 
ter diverges and had to be fixed manually, and the result 
is significantly worse than the FM curve) . As the Li curve 
has not been used extensively on local sightlines, it is dif- 
ficult to interpret the results, and we will not discuss it 
further. 

The results from our various fits are presented in Ta- 
ble [HI Note that despite the qualitative similarity of 
the curve to MW and LMC sightlines, three major pa- 
rameters {Rv = 4.17 ± 0.15, C3 = 1.70 ± 0.29, and 
C4 = 0.28 ± 0.07) differ significantly from the average 
values (7 is consistent with the average value). In gen- 
eral, the GRB 080607 sightline is UV-greyer, and its 
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Fig. 5. — Different extinction curves shown against our afterglow 
data, shown as the selective extinction E{X — V) = A\ — Ay- All 
curves are normalized to match the observed B — V color (tradi- 
tionally, UV extinction curves are plotted as E(\ — V)/E(B — V)). 
This illustrates the flatter nature of the derived extinction curve 
(higher Rv) and weaker 2175 A bump required along the GRB 
sightline relative to the average MW or LMC sightlines. 

2175 A bump weaker, than the average MW sighthnc. 
Still, all these properties are in the range seen along dif- 
feren t sightlines locally (e.g., from iFitzpatrick fc Massal 
[19901 : 2.3 < Ry < 6.6, 1.2 < C3 < 4.5, 0.15 < C4 < 0.90). 
No single local analog appears to match the properties 
seen toward the GRB exactly, but it is nevertheless no- 
table that our data are so well fit by the standard, locally 
derived laws without the need for any unusual parame- 
ters.^^ We will further discuss the implications of the 
FM parameters in §4.31 

Assuming any particular extinction model is not 
strictly necessary for this GRB: the large extinction 
column actually allows us to directly measure the 
wavelength-dependent extinction without need for fit- 
ting. Traditionally, UV extinction curves are presented 
as E{X — V) (i.e., A\ — Ay; the optical extinction itself 
Ay need not be known). Our if^-band measurement cor- 
responds to the rest-frame V band, and so if the intrinsic 
slope can be assumed, one can simply measure this value 
for each filter (or wavelength bin) by comparing the ob- 
served \-V color to the predicted color for the assumed 
intrinsic spectrum. The results are plotted in Figure [3 
illustrating the intrinsic differences between the curves 
and the inability of most of them to fit the data. 

3.4. Effect of Varying Intrinsic (3 

The above quoted results all assume (3 — 0.7. In reality, 
we do not know the exact intrinsic spectral index, which 
varies from burst to burst. As previously mentioned, 
the extinction of this burst is sufficiently large, and the 
intrinsic variation in /3 between events relatively small, 
that the errors introduced from variation in the spectral 
index are small. Here we quantify that statement and 
propagate the effects into our parameter uncertainties. 

This is not simply a matter of the flexibility of the fitting 
function: the model is quite limited in scope, with only four free 
parameters, each of which is constrained to a small allowable range. 
Indeed, some reported extragalactic sightlines, e.g., the high-2 
QSO sightline of ltvlaiolino et al.i (i2004( ). cannot be accurately fitted 
within this model. 



iKann et al.l (|2010f) have compiled photometry for a 
large number of bright, well-observed, Swift-era, GRBs 
and performed fits to the extinction (using the stan- 
dard MW/LMC/SMC method) and spectral index of 
each event. We downloaded the data in Table 2 of that 
work and removed all events which did not have a best- 
fit (among the three models) Ay < 0.2 mag within 2a 
to exclude events with significant or poorly determined 
extinction. We further removed any events reporting 
an unphysical Ay < mag at more than 2a and any 
event with an uncertainty in its derived spectral index 
ap > 0.2. The intrinsic spectral indices of this final sam- 
ple of 21 low-extinction, well-constrained bursts have an 
average spectral index of /? — 0.70 and standard devia- 
tion (T^ = 0.26. We take this as a representative sample 
with which to determine a prior on the intrinsic (unex- 
tinguished) spectral index /?. 

The observed spectral index between the J and 
bands for this GRB is (3 ~ 3.5, so the impact of red- 
dening (between these wavelengths) from dust is clearly 
much larger (by about an order of magnitude) than the 
typical variation in the intrinsic spectral index. This vari- 
ation in the intrinsic index is, however, the largest source 
of uncertainty in the measurement of the extinction pa- 
rameters. To take this into account, we refit our preferred 
extinction models for the ilcr cases and combined the re- 
sulting variation of the best-fit value in quadrature with 
the statistical uncertainties on the (3 = 0.7 fit. The final 
values for all extinction parameters (using the Reichart 
C2 — Ry correlation; the Fitzpatrick correlation is not 
significantly different) are presented in Table [51 

As an alternative to assuming an intrinsic optical /?, 
we also attempted our fits by including the X-ray flux 
value at the extraction epoch and assuming an unbro- 
ken power law over the full range between the optical 
and X-ray data (which allows for a much more precise 
derivation of /? as well as a constraint on the overall flux 
normalization, though it is strongly dependent on this 
assumption of an unbroken intrinsic index). This gives 
generally quite consistent values with our optical-only fit, 
in further support of our assertion that the derived dust 
properties are not strongly affected by our assumptions 
about the intrinsic spectrum. 

3.5. X-ray Scattering? 

Of particular not e in Table [51 and co nsistent with 
our previous work (jProchaska et al.l |2009[ ) , is the con- 
clusion of a large extinction column {Ay = 3.26 ± 0.35 
mag). This identification of GRB 080607 as a highly ex- 
tinguished event makes it a potentially useful test case 
o f the X-ray scat tering model for early-time afterglows 
((Shen et al.ll2009( ). However, even Ay = 3 mag is gen- 
erally inadequate to expect any significant effects on the 
X- ray light curve in this case. Following the discussion 
by iShen et all (|2009D . we calculate the 1 keV specific 
fluence f rom the prompt emission using the parameters 
given in iGolenetskii et al] (|2008[ ) and integrate the X- 
ray afterglow flux (starting at 100 s, and ignoring the X- 
ray flare) using our power-law fit. The resulting ratio of 
5'AG/S'prompt = (1-0 X 10-6)/(1.6 X 10"^) = 6.3 placcs an 
upper limit on the scattering opacity at this wavelength 
(i.e., Tscat < 6.3). Translating thi s to a limit o n the opti- 
cal opacity using Equation (8) of IShen et all ([2009) , the 
limiting dust extinction for this case is the thoroughly 
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TABLE 5 
Extinction Fits 



Model 


Av 


Ry 


Cl 


C2 


C3 


C4 


7 


xVdof 




(mag) 

















Average MW 

Average LMC 

LMC2 

SMC 

CCM 

FM+tie 

FM+Reichart 

Li 



1.25 ±0.03 
1.09 ±0.02 
0.16 ±0.03 
0±0.01 
0.82 ±0.06 

3.26 ±0.31 
3.52 ±0.35 
1.70 ±0.06 



3.1 
3.2 
2.6 
2.73 

2.41 ±0.12 
4.17 ±0.15 
4.69 ±0.19 



-0.07 
-1.28 
-2.16 
-4.96 

1.11 ±0.12 
1.29 ±0.15 
200 



0.70 
1.11 
1.31 
2.26 

0.31 ±0.04 
0.30 ±0.05 
12.3 ±0.6 



3.23 
2.73 
1.92 
0.37 

1.70 ±0.29 
1.66 ±0.30 
14 ± 285 



0.41 
0.64 
0.42 
0.46 

0.28 ±0.07 
0.31 ±0.07 
0.03 ±0.01 



0.99 
0.91 
1.05 
0.99 

1.10 ±0.06 
1.07 ±0.07 



127 / 24 
275 / 24 
1143 / 24 
1159 / 24 
123 / 22 
24.2 / 20 
22.9 / 20 
38.7 / 20 



Note. — Comparison of fits to the SED of GRB 080607 using a variety of extinction models, most of which 
cannot adequately fit the observations. Because the optical spectrum and photometry dominate the observations, 
most models converge to a low extinction value to try to accommodate the weak 2175 A bump and seemingly flat 
spectrum. These models are not consistent with the red IR color. Both a high Ry and a low C3 are required to 
explain the optical and IR data together, as reflected in the FM fits. Parameter uncertainties do not include the 
effect of the uncertain intrinsic spectral index P (a value of 0.7 is assumed). 



TABLE 6 

FM Extinction Parameters for 
GRB 080607 





Optical alone 


Optical + X-ray 


Parameter 


value 


value 


P 


0.7 ± 0.26 


1.08 ± 0.05 


Av 


3.26 ± 0.35 


3.07 ± 0.32 


Ry 


4.17 ± 0.25 


4.52 ± 0.23 


Cl 


1.11 ± 0.20 


1.37 ± 0.15 


C2 


0.31 ± 0.07 


0.22 ± 0.05 


C3 


1.70 ± 0.30 


1.82 ± 0.32 


C4 


0.28 ± 0.08 


0.37 ± 0.08 


7 


1.10 ± 0.07 


1.07 ± 0.06 


Xo 


4.596 





Note. — Final FM extinction parameters 
for GRB 080607. The values in the left col- 
umn incorporate only the optical data and 
include the effect of unknown intrinsic spec- 
tral index. Values at right assume an unbro- 
ken power law between the optical and X- 
rays. ci and C2 are tied to Rv as described in 
the text; Ry is significantly higher than the 
average MW or LMC curve but has a typical 
value for dense sightlines. The 2175 A bump 
(strength given by C3), ubiquitous in the MW 
but nearly absent in the SMC, is present but 
weaker than in the MW or LMC. 
unconstraining Ay < 686 mag, a value about 200 times 
higher than our direct measurement. Equivalently, the 
total fluence of the dust-scattered X-rays for this event 
is anticipated to be 200 times lower than the actual af- 
terglow fluence observed, and therefore undetectable. In- 
deed, for this event, the X-ray light curve follows a simple 
unbroken power law and (with the exception of th e early 
X-ray flare, a prompt-emis sion feature; iKocevski et al.l 
l2007t IChincarini et al.ll2007[ ) no significant h a rdnes s vari- 
ations of the type predicted bv IShen et alj ()2009t ). We 
conclude that, despite the large surrounding dust col- 
umn. X-ray scattering is not significant for this GRB. 

4. DISCUSSION 

4.1. Afterglow Luminosity in Context 

The impressive optical brightness for an event at z = 3 
has already been noted. In fact, as we shall show, 



after the effects of extinction are taken into account, 
GRB 080607 is among the most optically luminous GRBs 
(and therefore objects of any sort) to date, second only 
to GRB 080 319B. 

Following [Kann et all ()2007l |2010[) . we select as our 
comparison filter the z = 1 R band (that is, the wave- 
length which is shifted to the observed R band if at 
z = 1); this corresponds roughly to the rest-frame U 
band. For GRB 080607, this is shifted all the way to 
approximately the observed J band. Therefore, taking 
advantage of the apparent lack of color change, we shift 
all other filters to the J-band light curve using our model 
fluxes, extending this curve back to the observed emis- 
sion peak. This curve is corrected for Galactic extinction 
(only 0.02 mag), for host extinction (4.94 mag), and for 
the difference in luminosity distance between z = 1 and 
z = 3.036 using standard cosmology {h=0.7, ^Im = 0.3, 
JIa = 0.7). A small K-correction is then applied to match 
the spectrum exactly with the z = I R band, and the 
light curve is scaled (undilated) to z = 1. 

The result is plotted in Figure IH compared with the 
light curves of GRB 080319B and the three next mos t 
luminous events (from Figure 7 of iBloom et all 120091) . 
and with the peak luminosities of a la rge sample of wel l- 
studied Swift bursts (from Figure 7 of lKann et afllMol 'l 
At the beginning of observations GRB 080607 is compa- 
rable in luminosity to GRB 080319B (and at early times 
was likely brighter), but the prompt optical flaring of 
GRB 0803 19B pushes that burst to a higher luminos- 
ity over the next several minutes. GRB 080607 remains 
among the five most luminous bursts for the rest of its 
observed evolution. This illustrates the remarkable at- 
tributes of this burst that allowed it to provide such 
a detailed analysis of its environment. In terms of the 
afterglow (external shock) emission alone, GRB 080607 
may yet be the most luminous: GRB 080319B's peak 
appears to correlate with its prompt emission and fades 
particularly rapidly when the prompt emission ends; 
the ori gin of its early-time optical emission is still de- 
bated (jRacusin et al.l 12008'). The optical light curve of 
GRB 080607 bears no relation to the prompt emission 
and is certainly external shock-dominated at all times. 

4.2. Physical Properties 
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Fig. 6.— Near-UV luminosity of GRB 080607 (bold curve) com- 
pared to several other prominent bursts, as well as to a large sample 
of rapidly observed Swift GRBs from Kann ot al. (2010). Colored 
points indicate peak observed luminosities of the events described 
in that paper (unfilled points are events caught after the peak 
and therefore only lower limits on the peak luminosity.) At peak, 
GRB 080607 is among the most luminous GRBs known, peaking 
at Mu ^ —37 mag. At z = 1 it would peak at mag 6 if unobscured 
by dust. 
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Fig. 7.— X-ray luminosity of GRB 080607 (bold red curve) 
compared to all other 5™/4-followed GRBs. Several other promi- 
nent bursts are also individually colored. GRB 080607 is among 
the most X-ray luminous bursts at peak, but fades quickly to an 
average luminosity by later times. 



Unfortunately, the physical properties responsible for 
making GRB 080607 so energetic remain mostly hidden 
from view. In the burst rest frame, our optical obser- 
vations extend only to 10^ s and the X-ray observations 
cease at t « 1 day, which does not usefully constrain the 
jet opening angle. Conservatively setting tjct > 6 x 10"* s 
using the X-ray light curv e, following th e standar d equa- 
tions for the jetting time (|Sari fc Piran|[l999i : .Frail et alj 
1200 It) and fiducial values of density n = 100 cm" 
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This is an unusually large value of n, motivated by the ap- 
parent low value of the cooling break Uc and inference of a dense 
molecular environment along the line of sight, as discussed later 
in this section. Fortunately, the value of n only weakly alfects the 
derived value of 6 and E-.,. 



and efficiency rj = 0.2, we measure a jet opening angle 
^100 degrees; the equivalent lower limit 
on the beaming-corrected gamma-ray energy release is 
> 1.8 X lO^^f^Q^^^n^QQ erg, a fairly typical value. It is 
therefore not clear whether the extreme apparent lumi- 
nosity of this burst is attributable to intrinsically large 
energetics ( Cenko et al. 2010), favorable viewing angle 
(of a nonuniform, centrally c oncentrated jet, as w as sug- 
gested for GRB 080319B bv lRacusin et al.ll2008D . an in- 
trinsically narrowly concentrated (uniform) jet, or some 
combination of these parameters. 

To a large extent, the optical luminosity is simply an- 
other reflection of the total energe tics of the burst it- 
self: both in theorv (ISari et al.lll998l) and observationally 
(iGehrels et al.ll2008HNvsewander et al.ll2009HKann et al.l 
|2010D . the inferred afterglow luminosity scales approxi- 
mately linearly with E^iso, and if the optical light curve is 
extrapolated to late times the predicted optical flux is in 
the middle of the fluence-normalized distribution. How- 
ever, there is more to the story: the X-ray light curve of 
this burst is (except at the earliest times) not particularly 
bright; when normalized to the burst fluence it is quite 
typical for a Swift burst at early times (and actually is 
unusually faint at late times, due to its rapid unbroken 
decay) . 

There are two broad ways to interpret this. The sim- 
plest interpretation is that the cooling-break frequency 
Vc has a particularly low value compared to most GRBs, 
perhaps even below the optical band (the available data 
are marginally consistent with the X-ray and optical 
bands being on a single spectral power law). The obvi- 
ous culprit for this involves the external density n: The 
X-ray flux should be independent of density (assuming 
that vx > Vc for most bursts), but the cooling break and 
optical flux are sensitive to it [vc oc n^^; below the cool- 
ing break F oc n^/^). An external density 10 or 100 times 
the "typical" Swift value would push the cooling break 
from its typical position between the optical and X-ray 
bands into or below the optical band at early times, in- 
creasing the optical luminosity. Indeed, after correcting 
for extinction the early-time broadband SED appears to 
demand a low-cooling break: the optical-to-X-ray index 
at only 300 s is f3ox — 1-1, consistent with the X-ray 
spectral slope (/3x = 1.16 ±0.13). The probable low 
value of the cooling break also helps explain why a simi- 
lar extinction column is derived whether the optical data 
is considered alone (the most general case) or in conjunc- 
tion with the X-ray data assuming an unbroken power- 
law (which requires Vc < Vopt), as demonstrated in §3.41 
Unfortunately, the period of simultaneous temporal cov- 
erage between the optical and X-ray observations is too 
short to determine, via the light curve, whether a break 
is present between the bands. (The burst could also have 
exploded into a wind-stratified medium — one with vari- 
able density n oc — in which case Vc rises with time 
and the optical flux fades more rapidly than the X-ray 
flux, as is observed.) 

The alternative interpretation is that the optical 
flux originates from a separate emission component, 
the most obvious cand idate being the reverse shock 
eszaros fc ReesI |1997[ ). Several previous early fast- 
fading light c urves have been as s ociated with reverse 
shocks (e.g., lAkerlof et al.l 1 19991 : iKobavashi fc Zhangl 
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l2003tlPerlev et al.ll2008l: [Steele et al.ll2009D : qualitatively, 
the behavior of GRB 080607 appears similar to these 
events, although the decay is somewhat slower and there 
are no late-time observations to determine whether the 
light curve became forward-shock dominated as pre- 
dicted. The f actors determining the luminosity of the 
reverse shock (jZhang et al.ll2003f ) are generally the same 
as for the forward shock (and so a high external density 
would similarly aid the production of a luminous after- 
glow) , but can be further amplified if the magnetization 
Rb = £B,r/iB,f of the reverse shock is hi gh, due (for in- 
stanc e) to primordial fields in the ejecta (jGomboc et alJ 
[2008h . 

Unfortunately, the lack of late-time observations (to 
search for the appearance of a forward shock) or radio 
data (to more directly constrain n) prevents us from dis- 
tinguishing between these possibilities. Fortunately, we 
can speak more confidently about the other aspect of this 
burst's remarkable luminosity: the fact that it peaked so 
early (even if two bursts have similar energetics and late- 
time luminosities, the power-law nature of GRB light 
curves ensures that the event with the earlier peak time 
will have significantly larger peak brightness, fleeting as 
it is.) 

The peak time (for v > Vm) is set by the deceleration 
timcscale of the ejecta (jSari fc Pi ran 1999). Because the 
afterglow has already peaked and is fading at the start 
of our observations, the ejecta must have accumulated 
enough circumstellar matter to begin to decelerate and 
develop an external shock by this time: a mere 24 s after 
the BAT trigger (6 s in the rest frame), or more conser- 
vatively 32 s after the beginning of the prompt emission 
(8 s in the rest frame). 

Such rapid deceleration generally requires a high initial 
Lorentz factor F, although a very high interstella r den- 
sity also contributes. Using Equation 3 of Rvkof f~t al.l 

(|2009f )^^. we estimate F > 660?7g^2^nj^Qj/^, where 770.2 and 
nioo indicate values of the efficiency and external density 
relative to fidu cial values of 0.2 an d 100 cm~'^, respec- 
tively (see also iMolinari et aI]|2007D . Based on the pre- 
ceding discussion of the late-time optical luminosity, we 
have chosen an unusually large value for the interstellar 
density; even in this case the constraint on F is at the 
top end of the afterglow-inferred range (if still somewhat 
below the p air-opacity limits recently provided by the 
Fermi-hAT: lAbdo et al.ll2009D . It is notable that both 
Eiso and F are exceptionally large for this burst, which 
could suggest that the properties may be correlated. 

4.3. X-Ray and Optical Properties: the Environment of 
GRB 080607 

The derivation of precise values for the extinction pa- 
rameters along the GRB 080607 sightline (Table [5]) gives 
us an additional means for learning about its host envi- 
ronment. Although the reason for the variation of these 
parameters is not well understood even within the MW, 
some broad conclusions can be drawn. 

^'^ This equation is strictly valid only for the thin-shell sce- 
nario, in_whichthc burst duration is less than the deceleration 
time llSari fc Piran 1999; Mcszaros 2006). This is not strictly true 
for this GRB, as prompt emission is observed to continue during 
the light-curve decline. Howrever, as noted previously, the energet- 
ics are dominated by a single, bright pulse which ends well before 
the start of optical observations. 



First, we note the high value of Ry ~ 4 (or equiva- 
lently, since the parameters are tied in our modeling, the 
small value of C2), indicating a relatively flat extinction 
curve. In the diffuse interstellar medium (in the MW and 
in other galaxies as well), Ry typically takes on lower 
values of 2-4. UV-flat extinction curves are generally re- 
stricted to denser sightlines, proba bly because grains ar e 
able to coagulate to larger sizes ()Valencic et al.l [200I) . 
(At the same time, however, dense regions can also have 
low values of Ry as well as high ones.) The high Ry 
value is therefore suggestive of a dense environment — 
fully consistent with the conclusion from the atomic and 
molecular analysis that the sightline penetrates through 
a dark molecular cloud in its host. 

Second, the value of C3 is nonzero, indicating a signifi- 
cant 2175 A absorption bump. This is one of only a few 
clear detections of this feature at cos molog ical distances 
(iMotta et al.ll2002l:lJtinkkarinen et"aIll2004l:lEllison et al.l 
I2OO6I : iKriihler et all 120081 : lElfasdottiT et al .1 120091 ) and 
the highest-redshift detection of the feature yet. The 
identity of the carrier is still unknown (although poly- 
cyclic aromatic hydrocarbons and graphit e are consid- 
ered promising candidates; see lDrainel[2003l for a review 
and the processes that cause it to be present or absent 
are similarly not yet certain: an evolved stellar popu- 
lation (Noll et al. 2007), metalhcity (Fitz patrick 20M), 
the st rength of the UV radiation field (e.g.. lGordon et al] 
Il997f). and disturba nce of the environment due to shocks 
( Seab fc Shul]|ll983[ ) have all been cited in explaining its 
absence. Generically, however, it seems to be present in 
almost all sightlines in the MW and LMC, and in nearby 
disk galaxies, but absent in more disturbe d locations suc h 
as the SMC, nearby starburst galaxies (|Gordonll2005| ). 
and a t least one highly d isturbed sightline within the 
MW (jValencic et al.l l2003D . This suggests that the in- 
terstellar medium of the host of GRB 080607 is a closer 
analog of the more quiescent environments found in the 
MW and LMC than of the extreme conditions of nearby 
galaxies having high specific star-formation rates. Our 
upcoming study of the host galaxy (Chen et al. 2011) 
may help test this hypothesis. 

The strength of the bump is, however, weaker than in 
almost any sightline in either the MW or LMC (Figure 
[8|). Furthermore, the degree to which the bump is weaker 
does not follow the local correlations: in the MW, very 
low values of C3 tend to correlate with very low values 
of C2 (weakly) and 7 (strongly). In our case, a low C2 is 
observed, but it is still much higher than for the Orion 
Nebula sightlines in which the lowest values of C3 are 
seen. This may be an indicator that a different phe- 
nomenon is suppressing this carrier than is in operation 
within the MW Galaxy. Metallicity is not likely the cul- 
prit: the molecular cloud giving rise to the observed ex- 
tinction has near-Solar metallicity despite being at z > 3 
(jProchaska et al.ll2009t ). 

The strength of the FUV rise, C4, is fairly typical for 
local sightlines. However, the origin of the rise is even 
less secure than that of the 2175 A bump, and does not 
significantly constrain the environment. 

The X-ray inferred host-galaxy equivalent hydrogen 
column of GRB 080607 was measured to be A^h = 
cm , and comparable to the neutral hy- 
drogen column Nhi = 1.5^q'5 x 10^^ derived from the 
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Fig. 8. — Extin ction parameters for GRB 080607, compared to various MW and LMC sightlines from 'Va lencic et al.l 1*2004') and 
IMisselt et all l fl999l) . Diffuse MW sigfitlines are indicated with small, gray open circles; dense sightlines are indicated with small, filled 
black circles. Peculiar MW sightlines incompatible with the standard CCM one-parameter family are identified as blue circles. LMC 
sightlines are indicated with rectangles; the SMC curve is plotted as a star. Extinction parameter C2 is a measure of the UV slope (inverse 
greyness); C3 is a measure of the strength of the 2f75 A bump and 7 is a measure of its width; and C4 indicates the strength of the far-UV 
rise. All parameter values (and all pairs of two values) are within the distribution seen locally, although there is no single example of a 
local sightline that is consistent with the extinction properties of the GRB 080607 sightline in all aspects. 



damped Lyman-a line (|Prochaska et al.ll2009l ). This is 
a very large value, even considering the high extinc- 
tion in this direction: the ratio of Nii/Ay = 8 x 10^^ 
cm~^/mag is several times larger than observed in the 
MW, although quite typical of GRBs fo r which both 
value s have been securely measured (e.g., ISchadv et al.l 
120101 ). It is possible that this arises for reasons unre- 
lated to the molecular cloud — for example, if additional 
dust-free gas is located closer or further along the sight- 
line relative to the molecular cloud that is responsible for 
the absorption. If intrinsic, this combination of a weak 
(but present) bump and a large Nn/A y ratio is consis- 
tent w ith t he correlation discussed by lEhasdottir et al.l 
((2OO9I) and lGordoneFall (l200a) . 

5. CONCLUSIONS 

One of the brightest and best-studied GRBs (at early 
times) of the Swift era, GRB 080607 holds particular 
potential for revealing the nature of GRBs and their en- 
vironments at high redshift. While the relatively limited 
observed temporal range restricts our ability to study the 
intrinsic nature of this event, this is more than compen- 
sated by the abundant early-time optical/IR data that 
reveal the detailed properties of the dark-cloud sightline 
in its distant host. 

The utility of this event is perhaps most evident in the 
context of the class of "dark" GRBs. Many factors, both 
intrinsic (high i?iso, T, and n) and extrinsic (large but 
not extremely large Ay, a redshift placing the 2175 A 
bump in the optical window, and the fortuitous ability 
to observe immediately with telescopes in both the conti- 
nental US and Hawaii) had to conspire together to allow 
an event to be observable in such rich detail. Had this 
event been slightly less luminous ("only" comparable to 
GRB 990123, ^2 mag fainter at most epochs), its after- 
glow would have been only marginally detected, and only 
at the earliest times; further decrease in luminosity would 
have rendered it undetectable to small telescopes. Even 
a modest increase in the amount of extinction (higher by 
Av ~ 1-2 mag) or the presence of relatively UV-opaque 
SMC-like dust would have a similar impact, suppressing 
all of the optical measurements. 

The literature contains many examples of such sources: 



GRBs with a comparable dust column but insufficient 
luminosity to shine th rough it. Some prom inent cases 
include GR Bs 970828 ()Diorgovski et al. 2001). 06 0923A 
(iTaiivir et al...2008bl). 06 1222A an d 070521 (Cenko eF^ 
120091: iPerlev et al.ll2009D. 70306 (iJaunsen et al.ll2008D 
081221 ()Tanvir et al.ll2008aD . and 090709A ("C enko et al.l 
[20H . But even these objects were unusually bright or 
had particularly rapid or deep observations in their fa- 
vor. A truly typical-luminosity Swift event without rapid 
or deep observations would completely escape notice in 
most cases, permitting only shallow limits on its extinc- 
tion column. Therefore, there is every reason to think 
that very dusty environments like that of GRB 080607 
are actually not uncommon among GRBs (if not neces- 
sarily ubiquitous). This is in agreement with our after- 
glow plus host survey with the P60 an d Keck telescopes 
(ICenko et an[2009HPCTlev et al.l[2009t ). 

The extinction curve along the GRB 080607 sightline 
— a dark molecular cloud at z — 3.036 — is quite 
similar to that of our own Galaxy (with a significant 
2175 A bump), though there are differences in finer de- 
tails. The success in modeling the extinction curve of 
this event within entirely locally developed models is in 
some ways heartening, giving us confidence that with 
sufficient knowledge we should be able to understand 
the absorption properties even out to these immense dis- 
tances. As perhaps the most detailed determination of 
the extinction properties of a galaxy at cosmological red- 
shift to date, we suggest that the extinction curve in this 
work may be of use to others attempting to take into 
account the effects in other galaxies at high redshift (see 
Appendix). At the same time, some other GRBs and 
other techniques have also at times pointed to extinction 
curves that diverge dramatically from local templates, so 
the topic should continue to be addressed with caution. 

Once corrected for extinction, GRB 080607 rivals the 
"naked eye burst" GRB 080319B as the most luminous 
known object in the Universe. This extreme early lumi- 
nosity of GRB 080607 is likely the product of a variety of 
factors: it has one of the largest Eiso values to date, and 
its optical luminosity may have been further amplified 
by a large circumburst density in its host (or, alterna- 
tively, a bright reverse shock). Unlike GRB 080319B, 



Ultraluminous, Dusty GRB 080607 



13 



GRB 080607 has a smooth optical peak and shows no 
correlation with prompt emission at that time. Events 
like GRB 080607 demonstrate the power of GRBs to il- 
luminate the darkest corners of the Universe: not just 
the reionization era (on which much current attention is 
focused) but also the dustiest regions over the following 
several billion years when the global star formation rate 

— much of it occurring behind optically thick dust clouds 

— was at its maximum. Such dust-obscured regions are 
extremely difficult to study by other techniques, or even 
with most GRBs, as demonstrated by the class of "dark" 
bursts. The combination of early observations and ex- 
treme energetics of GRB 080607 were enough to over- 
come even this difficulty, and demonstrate the power of 
rare, individual events to illuminate these hard-to-study 
regions and improve our understanding of the early uni- 
verse. 
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APPENDIX 

APPLYING THE GRB 080607 EXTINCTION LAW 

Because the extinction law we derive for GRB 080607 is contained within the broader Fitzpatrick parameterization, 
application of this extinction law is relatively straightforward. Using the function f m.unred . pro in the GSFC IDL 
library, the extinction law can be calculated using, for example, the following simple code. 

function extinction080607, wav, rv, av 

; Returns the extinction (in magnitudes) as a function of rest wavelength 
; (in Angstroms) for the GRB080607 sightline, normalized to A_V. 

if n_elements (rv) eq then rv = 4.17 
if n_elements (av) eq then av = 1.0 
ebv = av/rv 

influx = replicated., n_elements (wav) ) 

fm_unred, wav, influx, -ebv, extflux, r_v=rv, c3=1.70, c4=0.28, gamma=1.10 
return, -2. 5*alogl0 (extflux) 

end 

Here Ry is a free parameter, left to vary in principle — although use of the default value is strongly recommended, 
since in this implementation it is tied to the Ci and C2 parameters within fm_unred.pro. 



